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ABSTRACT 

Existence  and  uniqueness  of  weak  solutions  are  shown  for  different  models  of  the 
dynamic  behavior  of  elastomers.  The  models  are  based  on  a  nonlinear  stress-strain 
relationship  (satisfying  a  locally  Lipschitz  and  affine  domination  property)  and  in¬ 
corporate  hysteretic  effects  as  well.  The  results  provide  alternatives  to  previous 
theories  that  required  monotonicity  assumptions  on  the  nonlinearities.  Results  with 
a  nonlinear  constitutive  law  and  nonlinear  internal  dynamics  are  presented  for  the 
first  time. 


1  Introduction 

In  this  paper  we  examine  the  theoretical  foundations  of  a  series  of  models  for  the  dynamic 
behavior  of  elastomers  (filled  rubber-like  materials).  As  outlined  in  previous  papers  [6,  7],  the 
basic  model  describing  the  longitudinal  motion  of  a  viscoelastic  bar  with  the  upper  end  ( x  =  0) 
fixed  and  a  tip  mass  on  the  lower  end  (x  =  £)  is  given  by 
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pAcutt  ~  {Aca)x  =  0  for  0  <x<£  (1.1) 

Mutt(t,£ )  +  Aca(t,£)  =  Mg  +  f(t)  (1.2) 

u(*,0)  =  0  (1.3) 

w(0,.T)  =  t/o  (1-4) 

ut{01x)  =  u1,  (1.5) 

where  p  is  the  mass  density,  Ac  is  cross-sectional  area,  M  is  tip  mass,  g  is  gravitational  accel¬ 
eration,  /  is  external  force,  and  u  denotes  the  stress.  If  there  is  no  tip  mass  and  Ac  =  p  =  1, 
then  we  can  write  the  system  in  variational  formulation  as 

Utt  oT  -  F(t)  ill  V*  (1.6) 

k(0,  x )  =  u0  £  V  (1-7) 

nt(0,  x)  =  ui  G  if,  (1.8) 


where  we  let  F  =  ffj((0,f),  H  =  T2(0,  ()  and  F(t )  =  f(t)fi#(x).  The  crucial  modeling  question  is 
what  type  of  stress  strain  relationship  best  describes  the  material.  In  our  earlier  experimentally 
based  investigations  of  elastomers  we  found  that  a  nonlinear  constitutive  law  is  needed  ([4,  5]), 
i.e., 


u{ux)  =  Eux  +  ge{ux)  +  CDutx , 


(1.9) 


where  the  third  term  on  the  right  side  is  a  first  approximation  to  a  damping  that  these  materials 
exhibit.  (For  simplicity  of  presentation  here  we  will  later  make  the  assumption  that  E  =  1.) 
Comparing  the  actual  experimental  data  with  this  model,  we  demonstrated  good  agreement  for 
unfilled  and  lightly  filled  rubber  bars.  However,  the  model  is  not  adequate  to  describe  medium 
or  highly  filled  elastomers  that  exhibit  significant  hysteretic  behavior.  To  account  for  this 
property  of  the  material  we  included  a  Boltzmann  integral  term  in  the  constitutive  relationship 


^"(^a;)  —  Eux  T  ge{Ux)  T  C JjUtx  T  J  &  ^  ^  gvi^x^  utx)ds. 


(1.10) 


As  described  in  [5]  under  specific  assumptions  on  the  prehistory  of  the  motion  this  relationship 
can  be  expressed  in  an  equivalent  internal  variable  formulation 


a(ux) 

£lt  +  Ci£i 

£ i ( 0 ,  x) 


Eux  +  ge{ux)  +  C'DUtx  +  £i 

~gjdv  {.'M’x\ utx) 

0. 


(1.11) 

(1.12) 

(1.13) 


We  can  think  of  £i  as  an  internal  strain  variable  whose  dynamics  is  described  by  (1.1 2)- ( 1 . 1 3) . 
We  also  note  that  although  these  constitutive  relationships  are  expressed  in  terms  of  the  in¬ 
finitesimal  strain  uXl  the  formulation  is  equivalent  to  considering  finite  strains  £  =  ux  +  |u2, 
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with  different  nonlinear  functions  ge  and  gv .  The  internal  variable  system  ( 1 . 1 1 )-( 1 . 13)  can  be 
further  generalized  by  considering  nonlinear  internal  dynamics,  i.e. ,  (1.11)  with 

£i t  ^rhi  ,  fir  {w’xi  'Utx')  T  5fm(^-i)  (1-14) 

at 

£l(0,  ®)  =  0.  (1.15) 

We  note  that  the  resulting  constitutive  law  is  no  longer  equivalent  to  a  Boltzmann  integral 
formulation.  A  further  generalization  that  has  proved  important  for  highly  filled  elastomers 
involves  multiple  internal  variables.  That  is,  one  replaces  £\  in  (1.11)  by  a  finite  sum  J2£j  °f 
internal  variables  £j  which  satisfy  systems  of  the  form  ( 1 . 12)-(  1 . 13)  or  (1.14)-(1.15). 

In  this  note  we  consider  the  well-posedness  of  the  basic  model  with  the  three  different 
constitutive  laws.  Thus  we  first  show  that  the  nonlinear  problem  with  no  hysteresis  ( ( 1 .6)- ( 1.8) 
with  (1.9))  is  well  posed  under  rather  general  assumptions  on  the  nonlinear  function  ge.  Our 
first  result  guarantees  the  existence  of  a  unique  local  weak  solution  under  a  local  Lipschitz 
condition  on  the  nonlinear  function.  If  we  impose  an  additional  growth  assumption  on  the 
nonlinearity,  then  the  weak  solution  is  global  (i.e.,  exists  for  any  time  interval  [0,  T]).  We  then 
show  that  similar  results  can  be  obtained  for  the  nonlinear  problem  (i.e.,  nonlinear  gt- .  gv )  with 
linear  internal  strain  dynamics  ( ( 1 .6)-(  1 .8)  with  ( 1 . 1 1 )- ( 1 . 13))  and  for  the  nonlinear  problem 
with  nonlinear  internal  strain  dynamics  ( ( 1 . 6)- ( 1.8)  with  ( 1 . 1 1 ) , ( 1 . 14)- ( 1 . 1 5) ) .  The  first  two 
problems  were  previously  studied  under  certain  monotonicity  and  growth  assumptions  on  the 
nonlinearities  (ge  and  gv)  in  [2]  and  [6],  respectively.  The  results  in  this  paper  demonstrate 
that  well-posedness  can  be  achieved  under  relaxed  assumptions  and  can  be  extended  to  the 
problem  with  nonlinear  internal  dynamics.  The  techniques  we  use  were  successfully  employed 
to  establish  existence-uniqueness  of  weak  solutions  for  linear  evolution  equations  of  second 
order  in  t  in  [8]  and  for  semilinear  second  order  evolution  equations  where  the  nonlinear  forcing 
term  satisfies  a  global  Lipschitz  condition  in  [9].  In  [1]  such  techniques  were  extended  to  study 
a  nonlinear  beam  equation  where  the  nonlinearity  satisfies  only  a  local  Lipschitz  condition. 
Arguments  for  our  results  concerning  the  nonlinear  problem  with  no  hysteresis  are  very  similar 
to  the  ones  used  in  [1],  We  give  a  fairly  detailed  exposition  here  in  order  to  be  able  to  easily 
refer  to  these  as  we  extend  the  well-posedness  result  to  the  systems  with  internal  variables. 


2  The  nonlinear  problem  with  no  hysteresis 

In  this  section  we  investigate  the  system 

Uft  (-'D^txx  UXx  (ffe{^x^x  —  F(t)  m  I  (2.16) 

u(0,  &)  =  no  £  V  (2.17) 

v/<(0,  x)  =  ui  £  H,  (2.18) 

where  V  =  L7^(0,f!)  >  H  =  X2(0,/)  ►  V*,  Cd  >  0  and  E  has  been  normalized  only  for  the 

sake  of  convenience.  We  make  the  following  assumptions: 
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(Ag'J  The  nonlinear  function  ge  satisfies  the  following  local  Lipschitz  condition:  let  H,  (0)  denote 
the  ball  of  radius  r  centered  at  0  in  if  and  for  some  positive  constant  LBr  we  have 

\\9e(w)  ~  5e(^)||  <  LBr\\w  -  v|| 

for  all  w,  v  G  Br(  0). 

(24e)  There  exist  constants  C\,Ci  such  that 

||<7e(*f)||  <  Ci||u||  +  C2, 

for  every  v  G  II . 

( Af )  The  forcing  term  F  satisfies 

F  £  L2(0,T-,V*). 

We  clehne  the  space  of  weak  solutions  to  be 

W(0,  T)  =  {u  G  L2(0,  /’:  V)  I  ut  G  L2i 0,  T:  V),  utt  G  L2(0,  V:  V*)}, 
with  norm  given  by 


Definition  2.1  We  define  u  G  i/(0,T)  to  be  a  weak  solution  of  (2.16)-(2.18)  if  it  satisfies 

t*.’  F)  r  * ,  k  t  T  T  {B]^p')v*,v  (2.19) 

for  every  p>  G  V 

and 

u( 0,  x)  =  u0  G  y  (2.20) 

ut(0,  x)  —  ui  G  H.  (2.21) 

Here  (■,  ■)  denotes  the  inner  product  in  if  while  (•,  -)v*,v  represents  the  usual  duality  product. 
We  first  prove  the  following  local  existence  theorem: 

Theorem  2.1  Under  assumptions  (As&)  and  (Af)  and  for  any  u0  G  V,u i  G  if,  the  system 
(2.19)-(2.21)  has  a  unique  weak  solution  on  some  interval  [0,t*]. 

Proof:  Let  P  denote  the  Hilbert  space  radial  retraction  onto  the  ball  in  H  centered  at  u0x  with 
radius  1,  and  let 

//.  ('•)  =  9e(Pv). 
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Thus  ge  satisfies  a  global  Lipschitz  condition 

11/7.  -  &WII  <  ,+ll-c, I:)  II “  Prll  -  L\\w  -  V 

for  every  w,v  G  H  and  also 

||<7e(w)||  <  L\\w\\  +  C, 

where  G  =  1 1  <7e ( 0 )  1 1  >  0. 

Now  we  consider  the  problem 

(■  1)1!  ^ xx  (/G  (T'.t;)  )x  F(l)  m  h 

u(o,  x)  =  uq  g  y 

ut(0,  x)  =  w-|  G  77. 

Let  {'(/;,}  ■'h-i  be  any  linearly  independent  total  subset  of  V.  For  each  m  let 

V"‘  =  spanjV’i,  i> 2,  ■  ■  ■ ,  #m}, 


(2.22) 

(2.23) 


(2.24) 


(2.25) 

(2.26) 
(2.27) 


and  choose  G  Vm  such  that  i/q  — »  u0  in  V  and  vf  — >  u i  in  H  as  to  — >  oo.  We  develop 

the  standard  Galerkin  approximations  for  the  problem  (2.25)-(2.27).  Let  um(t )  =  C’"  (/•)</',• 

be  the  unique  solution  of  the  m-dimensional  ordinary  differential  equation  system 


(utt^j)v*tv  +  Cd{u™xi  tii)  +  +  (&(«”)* ipjx}  =  (F,  ipj)v*,v  (2.28) 

um(0,x)  =  u™  (2.29) 


i(0,^)  =  <. 


(2.30) 


To  obtain  an  a  priori  estimate  we  multiply  (2.28)  by  and  sum  over  j  to  arrive  at 

{n* ,«ri.-,v  +  CB(C»3  +  (OS)  +  «.«).» S>  =  (ZOv.v-  (2.31) 

Thus 

1  d 
2di. 


{lK(i)H2  +  |K(*)H2}  +  CdIKIP  =  -<ftW‘),«S}  +  (2.32) 

Integrating  from  0  to  t  we  obtain 


IKWII2  +  IK(0H2  +  Wn  /  I|^(r)||2dr  =  ||< 

Jo 

+2  J  (F,  u™)v*,vdT. 


m  1 1 2  ,  ||  m  1 1 2 


+ 


Oar  I 


;i’(W 


m  \  m  \  i 
ux  ),UTX)dT 


(2.33) 


Using  assumption  (/Ij).,  the  boundedness  property  (2.24)  and  standard  inequalities,  we  obtain 

rt  __  __  L2  ft 


im2  +  \Km2  +  v  l  iK(r)n2dr  <  iKir  +  ikji2  +  T  /  ikmii2^ 

JO  0  Jo 


+ 


1  U 


/  ||F(t)||^t  + 
Jo 


C2T 


(2.34) 
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where  6  is  chosen  such  that  v  =  26'/)  —  36  >  0.  By  applying  the  Gronwall  inequality  we 
can  conclude  that  the  sequence  {||K||2}  is  bounded.  Hence  there  exists  a  positive  constant 
C  —  C(ui,u0,  ge,  F,T)  independent  of  m  such  that 

IKWII"  +  ll<(<)ir  +  V  f  ll«T)||Jrfr  <  C.  (2.35) 

JO 

Now  we  can  argue  as  in  [1,  2,  3,  8]  that  there  exists  a  subsequence,  again  denoted  by  {a™},  and 
limit  functions  u  £  W1,2{ 0,  V:  V )  and  g  £  L2{ 0,  T;  77),  such  that 

um  — >  v/  weakly  in  H  1  2 ( 0.  V:  U)  (2.36) 

ge(u™)  g  weakly  in  X2(0,T;77).  (2.37) 

Using  these  convergences  we  obtain  that  u  satishes 


<««,  +  CD(utx,  <px)  +  («*,  (fx)  +  (g,  tpx)  =  (jp,  <f)v*,v  (2.38) 

for  every  9?  £  V  and 

•</(0,#)  =  w0  6  V  (2.39) 

v/<(0,  x)  =  tti  £  H.  (2.40) 

We  note  that  as  in  [1]  we  have  that 


is  continuous  over  P(0,  V:  V)  equipped  with  the  topology  of  L2(0,  T]V)  and  thus  by  density 
over  L2(0,T;V).  So  utt  £  (L2(0,T;V))*  =  L2(0,T;V*)  and  since  we  already  established  that 
u  £  TU1,2(0,  7 ;  V)  we  can  conclude  that  u  £  77(0,  T).  By  [8,  Remark  1,  p.  555]  this  also  provides 
the  additional  regularity: 


u  £  C([0,T];y)}.  uteC([0:T];H). 

To  conclude  that  u  is  indeed  a  weak  solution  of  (2.25)-(2.27)  we  need  to  show  that 

(g.  ifx)  -  (g,  (</,.),  95..;;)  for  every  9?  £  V.  (2-41) 

This  is  achieved  by  establishing  the  strong  convergence  u™(t)  — >  ux(t )  in  H  as  m  — >  00.  Let 
zm(t )  =  um(t)  —  u(t).  Using  u"'  and  ut  as  test  functions  in  (2.28)  and  (2.38),  respectively  and 
integrating,  we  find 

Ikrwil2  +  IKWf  +  2 CD  f  K(r)||2dr  =  IK  -  ^ll2  +  IK  -  tio-ll2 

JO 

-2/  <^e(«”(r))  -5(T)fi?™(T))^T  +2  [  {F(t),z™(t))v*,vcIt  +  (2.42) 

JO  JO 
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where 


Now 


xm(t)  =  2\-{ut(i),u?(t))  -  (ux(t),u%(t))  -2Cd  f {utx(t),  u™x{r))dT  +  (u1:u™) 

Jo 

+  <r'  [  (ge{u™(T)):uTX(T))dT  -  f  (g(r) ,  u™(r)). dr 

Jo  Jo 

+  2  J  {?’(-).  vt(t))<It  .  (2.43) 

2/  {9eOC{r))  -~g{T),z™x{T))dT  =  2  /  -  ^e(«e(r)), 

JO  JO 

+  2  jo(ge{ux(T)) -g(T),z?x(T))dT,  (2.44) 


and 


1 2  /  (//.  «(T) )  <1  {"At))..  z?x{t))<It\  <2  f  /.||<(t)  v/,r(-)||||-^.(-)||dr 

JO  JO 


< 


X  f  \K(r)\\2dr  +  6  f  lie (t)|| 2dr.  (2.45) 
0  Jo  jo 


Substituting  (2.45)  into  (2.42),  we  obtain 


kwh2  +  kwh2  +  (ac*  - «)  f  kooiimt  <  ik  -  «iii2  +  i«  -  ^n2 

Jo 

+T L  H'-MH2*'  +  NM<)I  +  in.(«)l.  (2.46) 

where 

Ym(t)  =  2  l  (ge(ux(T))  -  g{T),z™(r))dT  +  2  f  {F{r),z™{r))vydT. 

Jo  Jo 

Choosing  the  same  subsequence  as  before  we  have  ||it™  —  Wi||  — *  0,  || u™x  —  Wo®||  ~ ^  0.  Since 
z™  — »  0  weakly  in  X2(0,  T)  V )  we  also  have  that  |  Y„, ( /,)  |  — *■  0.  We  can  argue  that  |Wm(t)|  — ►  0  as 
m  — »  oo  for  a.e.  t  in  the  following  way:  using  the  convergences  (2.36)  and  (2.37)  in  (2.43)  we 
obtain  the  integrated  form  of  (2.38)  with  p>  =  ut.  The  fact  that  u  satisfies  this  equation  gives 
the  required  result.  Now  by  applying  the  generalized  Gronwall  inequality  to  (2.46)  we  can  see 
t  h  at 

lkrj)l|2  0  for  a.e.  t  e  [0,T], 

and  we  can  conclude  that 


-»■  )  strongly  in  h2(0,  7';  //). 

This  guarantees  (2.41)  and  thus  u  is  a  weak  solution  of  (2.25)-(2.27). 

Uniqueness  of  the  weak  solution  of  (2.25)-(2.27)  can  be  shown  in  the  standard  way  (e.g., 
see  [2,  3,  6,  8]). 
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We  now  prove  that  (2.16)-(2.18)  has  a  unique  weak  solution  on  some  interval  [0,f*].  By 
the  above  remarks  we  know  that  the  weak  solution  of  (2.25)-(2.27)  has  the  property  that  ux  is 
continuous  in  t.  Thus,  there  exists  t*  with  0  <  t*  <  T  such  that 

||ua.(t)  —  u0a.||  <  1  for  all  t  G  [0,U], 

and  therefore 

<1,  (".,  (/■))  =  ge{ux{t))  for  all  t  G  [0 ,t*]. 

Hence  u  is  a  weak  solution  of  (2. 16)-(2. 18)  on  [0,t*].  Uniqueness  of  the  weak  solution  can  again 
be  shown  in  the  standard  way.  This  completes  the  proof  of  Theorem  2.1. 

Now  we  use  the  additional  assumption  to  guarantee  the  existence  of  a  global  weak 

solution. 

Theorem  2.2  Under  assumptions  (A.)k),  (/W),  (Ar)  the  system  (2.16)-(2.18)  admits  a  unique 
global  weak  solution. 

Proof:  As  before  we  can  define  Galerkin  approximations  um(t)  =  C'_; (/)(/•;  to  solve  (2.28)- 

(2.30)  with  the  nonlinear  function  ge  instead  of  ge.  By  assumption  (A{,e)  we  can  develop  a  similar 
a  priori  bound: 

IKWir  +  \K(t)\\2  +  V  ['  \KA.r)\\2dT  <c=  C(u0,  u1:  F ,  T,  c1,  C2).  (2.47) 

Jo 

Thus  we  can  again  obtain  convergences  (2.36)-(2.37).  Additionally,  as  in  [2,  Lemma  5.1]  we 
can  show  that 

(/"’(/.)  — ►  ux(t)  weakly  in  H. 

(The  arguments  to  obtain  this  convergence  in  [2]  depend  only  on  the  a  priori  bound  and  the 
general  Arzela-Ascoli  Theorem  and  are  independent  of  the  specific  assumptions  on  the  nonlinear 
function.)  Thus  by  the  weak  lower  semicontinuity  of  the  norm  in  H  we  obtain  that 

kwh2  <  c. 

So  the  proof  can  be  completed  exactly  as  before  using  the  local  Lipschitz  property  of  ge  in  the 
ball  B 0)  in  H. 

3  The  nonlinear  problem  with  linear  internal  dynamics 


In  this  section  we  consider  the  system 

Uti  h'D^txx  ^ xx  (gei^x^x  &lx  B(t)  m  V  (3.48) 

.1  /  T  ©i€|  di  flv  (  ^ x :  U  !x )  (3.49) 

e^O)  =  0  (3.50) 

u(0,x)  =  u0  (E  V  (3.51) 

ut(0,x)  =  ui  G  H,  (3.52) 
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where 


fjvi  i^'x)  if  'U'tx 

flvd\  ^x  j  if  hla;  ^  0 


i.e.,  the  viscoelastic  response  function  is  different  when  the  strain  is  increasing  and  when  it  is 
decreasing.  The  internal  dynamics  is  interpreted  in  the  sense  that 


£i(f)  =  /  e  Cl[t  s)-^-gv(ux,usx)ds.  (3.53) 

Jo  as 

Integrating  (3.53)  by  parts  we  obtain 

Ci  e  gv(ux,usx)ds 

o 
K 

+  J2  h(f  ~  tk)e~Cl[t~tk\-l)k+1  [ gm{ux(tk ))  -  gvd{ux(tk))] ,  (3.54) 

k= 0 


where  h  is  the  Heaviside  function  and  tk ,  k  >  1,  are  the  points  where,  roughly  speaking, 
utx(tk )  =  0  with  t0  =  0  (e.g.,  see  [6]).  More  precisely,  as  explained  in  detail  in  [5],  [6]  the 
definition  of  gv  is  based  on  the  a  priori  given  set  of  points  {tk  {  where  the  value  of  gv  takes 
alternate  values  gm  and  gvd  on  successive  intervals  [tk:tk+ 1).  That  is,  in  the  system  formulation 
(based  on  experimental  data)  one  is  given  functions  gvi  and  gv#  and  a  sequence  of  points  {//,  ] 
so  that  gv  is  defined  by  the  alternating  values  gvi,gvd  on  intervals  [tk:tk+ 1).  Thus  in  essence  gv 
depends  on  ux  and  t  and  not  on  utx.  In  what  follows  we  will  therefore  use  the  notation  gv(ux). 


Definition  3.1  We  define  (u,e i)  G  U(0,  T)  x  T2(0,  T;  H)  to  be  a  weak  solution  of  (3-48)-(3.52) 
if  it  satisfies 


{utti  v}v*,v  +  CD(utx,  (px)  +  (ux,  <p>x)  +  ( ge{ux ),  <Px)  +  (ei,  <fx)  =  (F,  fi)v*,v  (3.55) 
for  every  G  V 

u(0,3i)  =  u0  G  V  (3.56) 

y/f (0.  ,r)  =  u\  G  H  (3.57) 


ft 

£i (t)  =  //,  (".,(/■))  -  /  gv(ux)ds 

J  o 
K 

+  Ht  ~  h)e~Cl(t~tk\-l)k+1  [ gvi{ux(tk ))  -  gvd{ux(tk ))] .  (3.58) 

k= 0 

We  make  similar  assumptions  on  as  on  ge,  namely, 

[Agf)  The  nonlinear  functions  gvi,gvd  satisfy  the  following  local  Lipschitz  condition:  let  Hr(0) 
denote  the  ball  of  radius  r  centered  at  0  in  H  and  for  some  positive  constants  L1Bt  and 
In  we  have 

JD  t 

|| gvi(w)  -  5«(«)||  <  LBt\\w  -  u|| 
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II 9vd(w)  -  9vd{v) ||  <  LdBr\\w  -  v|| 

for  all  w,  v  £  Br(  0). 

(Abv)  There  exist  constants  O^Oj^Oj  and  C %  such  that 

||5m(^)||  <  q||n||  +  Cj, 

||3W(^)||  <  C|||^|  +  C'j. 

for  every  v  £  II . 

Theorem  3.1  Under  assumptions  {A:llt).  (Agv),  ( A$ )  the  system  (3-48)-(3.52)  has  a  unique  local 
weak  solution. 


Proof:  The  proof  is  essentially  the  same  as  the  proof  of  Theorem  2.1,  so  we  just  outline  the 
crucial  steps.  First,  as  in  [6]  we  consider  the  interval  [0,  ij].  Let  P  denote  the  Hilbert  space 
radial  retraction  onto  the  ball  in  H  centered  at  u0x  with  radius  1,  and  define 


9e{v) 

=  ge{Pv ), 

(3.59) 

0vi(v) 

=  9vi{Pv ), 

(3.60) 

9vd{<') 

s> 

II 

(3.61) 

Thus  ge ,  gvi ,  gvd  satisfy  the  following  global  Lipschitz  and  boundedness  properties: 


h(w)  ~  &(w)||  <  Lj\\w  ~  «||,  3  =  G  vi-  vd\ 

ll&MII  <  ci\M+-ci 


(3.62) 

(3.63) 


Hence  we  can  consider  the  problem 

C  D^txx  ^ xx 

-  S-  gv{ux)  -  e~Cltgv(u0x)  -  Cl  f  e~Cl(t~s)gv{ux)ds  =  F(t)  in  V* , 
ox  L  Jo  J 

v/(0,  x)  =  no- 
ut{ 0,  x)  =  U) . 


(3.64) 

(3.65) 

(3.66) 


We  develop  Galerkin  approximations  {um}  to  (3.64)-(3.66).  The  additional  terms  in  this  system 
as  compared  to  (2.25)-(2.27)  cause  no  difficulties  in  obtaining  an  a  priori  estimate  similar  to 
(2.35)  due  to  the  properties  (3.63).  Thus  the  convergences  (2. 36)- (2. 37 )  can  be  obtained  and 
the  strong  convergence 

«”(<)  — >  ux(t )  in  H 

can  be  established.  One  can  also  easily  observe  the  additional  regularity:  u  £  0(0,  H;  V'"),  ut  £ 
( ’(0,  /. j ;  II ).  and  also  that  £i  £  0(0,  /p  //).  By  the  continuity  property  of  ux  in  t  there  exists 
I*  >  0  such  that  for  t  £  [0,0],  ||u;cr(^)  —  n0ar||  0  1.  Hence  gj(ux )  =  gj(ux),  j  =  e,vi,vd  on 
[0,0].  So  u  is  a  weak  solution  of  (3. 48)- (3. 52)  on  the  interval  [0,0].  If  O  =  tj  it  makes  sense 
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to  consider  the  next  interval  where  local  existence  of  a  weak  solution  can  be  established 

exactly  as  before.  Uniqueness  of  the  weak  solution  is  shown  in  the  standard  way. 

The  existence  of  a  global  weak  solution  can  be  guaranteed  under  additional  boundedness 
assumptions. 

Theorem  3.2  Under  assumptions  (Age),  (Agv),  (Ahe),  (Abv),  {Af)  the  system  (3.48)-(3.52)  ad¬ 
mits  a  unique  global  weak  solution. 

Proof:  First  we  consider  the  interval  [0,U]-  We  develop  Galerkin  approximations  {(/'"' }  for 
(3.48)-(3.52)  as  before  and  by  the  boundedness  properties  (vhe),  (zUJ  we  can  obtain  an  a  priori 
estimate  like  (2.47).  The  crucial  step  in  the  proof  is  that  we  can  establish  the  convergence 

u™(t)  — >  ux(t )  weakly  in  H, 

and  thus  guarantee  that  ||w.,.(/.)||2  <  C .  Now  the  local  Lipschitz  property  of  ge  and  gv  can  be  used 
in  the  ball  B 0)  to  yield  the  strong  convergence  u™(t)  — *  ux(t )  in  H  as  before.  Uniqueness 

can  be  shown  in  the  standard  way,  and  then  the  weak  solution  can  be  extended  to  the  next 
intervals  [U,L+i]  ,  i  >  U 


4  The  nonlinear  problem  with  nonlinear  internal  dynam¬ 


ics 

We  consider  the  system 

ytt  y^xx  (*7e  &ix  —  7^(t)  m  (4.67) 

£it  +  C i£i  =  gin(e i)  +  — gv(ux ,  utx )  (4.68) 

£i(0)  =  0  (4.69) 

u( 0,x)  =  uo  G  V  (4.70) 

i*t(0,  sc)  =  ui  €  H,  (4-71) 


where  the  nonlinear  functions  ge,gv  are  as  in  Section  3.  We  interpret  the  internal  dynamics  in 
the  sense  that  the  internal  strain  S\  solves 

Jrt 

e-ci(t-d  [-Clgv(ux)  +  gin{&ij]ds 
o 
K 

+  X)  h(t  -  tk)e~Cl(t~tk'>(-l)k+1  [gvi(ux(tk))  -  gvd(ux(tk))] ,  (4.72) 

k= 0 

where  again  h  is  the  Heaviside  function  and  U,  k  >  1  are  dehned  as  in  Section  3.  We  suppose 
that  the  nonlinear  functions  gK.  gv.  gvn  satisfy  global  Lipschitz  properties: 
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(Al)  For  some  positive  constants  Le ,  Lm ,  Lvd  and  Ltn  we  have 


II  9e(w) 

-Se(v)|| 

< 

Le||u;  —  v  || 

(4.73) 

\\9vt(W) 

-  9vi(v)  || 

< 

Lvi\\w  -  v|| 

(4.74) 

II  9vd(w) 

-  9vd(v)  || 

< 

Lvd\\w~v\\ 

(4.75) 

IlflwM 

-  ^n(»)|| 

< 

Lm\\w-v\\ 

(4.76) 

for  all  w,  v  £  H. 

We  prove  the  following  theorem: 

Theorem  4.1  Under  assumptions  (Al),  ( Ar)  the  system  (4-67)-(4-71)  has  a  unique  global  weak 
solution  . 


Proof:  Let  ns  first  consider  the  interval  [0,  £1].  On  this  interval  gv  =  gvi  or  gv  =  gv$  depending 
on  the  initial  conditions  and  the  forcing  term.  We  define  the  following  approximate  sequence 
{uN,e#}  !  let  be  tho  unique  weuk  solution  of  the  system 

iutt^)v*,v  +  Cniu^tpx)  +  (uZ,<p>x)  +  (ge(u^),<p>x)  +  (of.  ifx)  =  (/•'.  (4.77) 

uN( 0,.t)  =  u0  £  V  (4.78) 

u^(0,x)  —  ui  £  H  (4.79) 

and 

Ei  =  9v(ux)  -  e~Cltgv(u0x)  +  e_Cl(t_s)  [-d gv(ux)  +  ^„(ef _1)]  ds  (4.80) 

where  =  0.  The  sequence  is  well-defined  since  (4.77)-(4.80)  is  the  same  as  (3.55)-(3.58)  except 
for  the  known  term  /0  e~Cld~slgin(e^~  1(s)) ds.  It  is  easy  to  see  that  Theorem  3.2  extends  to  this 
case  as  well.  Also  the  assumption  (  Al)  for  ge  and  gv  guarantees  that  (Agi),  ( A:h ).  (/h,,)  and 
(,4ft,  )  are  satisfied.  Our  goal  is  to  show  that  {uA  (t)},  {ux(t )},  {eA (t)}  are  Cauchy  sequences 
(uniformly  in  t  <  tj). 

Let  uN  =  uN  —  uA_1  and  =  cA  —  eA_1.  Then  we  have 


,  <p)  +  CD(uZ,  <fx)  +  (uZ,  yx)  +  (ge(uZ) 


9e(uZ  1),¥x)  +  (ti  ,<Px)  =  o  (4.81) 


With  9?  =  uf  this  gives 

~{||if(«)||!  +  ||<(()in  +  cD||fi»||J  <  (L>  +  n)||«»(<)||||fiS(()|| 

+d 1*  f  ll*7«|l<fa|k"(<)ll  +  ('  ll^’M  -  ^-!(^)ll*l|iS(OII.  (4.82) 

Jo  Jo 
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where  Lv  =  max {//'*,  Thus, 

II^WII2  +  ||*?(*)||2  +  Wn  f  ||  ul(T)\\2dr  <  2(Z*>  +  r)  f  ||fif(T)||||«^(r)||dr 

JO  JO 

+  2cjLv  f  f  \\u* (s)||da||«^.(r)||dr 

JO  JO 

+  2 Lm  [  /  ||ef_1(-s) -ef“2(3)||^||h^(r)||Jr.  (4.83) 

JO  JO 

We  estimate  the  second  term  on  the  right  side  as 

2ci Lv  f  j  \\u* (s)\\ds\\u?x(T)\\dT  <  2eyLv  l  \\u* (r)||dT  f  ||«^(r)||dr 

JO  JO  JO  JO 


si  o 
c1Lvt1  /•' 


(r)||hr  +aLvS{ 


\Kx(T)\\dT 


jr1 1  ll^(T)H2^T  +clLvSlt1j0  ||«^(r)||2dr. 


(4.84) 


The  third  term  can  be  estimated  similarly  to  yield 


/  /  ll£f  1(-s)-£f  2{s)\\ds\\UrX{T)\\dT  <  ^-FT  /  Ikf  1(T)-£1  2(r)\\2dr 

JO  JO  0-^  Jo 


+  Vn8- 


%  r  n«£(i 

Jo 


(4.85) 


Hence  (4.83)  gives 


n«r(<)ii2  +  ii^(oir+/*i/  ii^(t)ii2^t  <  &(  ii^mii^t 

JO  JO 

+  M3  /t||£f-1(T)-£f-2(r)||2dr,  (4.86) 

Jo 

where  [ii  =  L  +L  1,1 ,  f/,3  =  L  s-2tl  and  hi  is  chosen  such  that  /i-\  =  26'/;  —  ( //  +  iT)h2  — 
Cli/h2H  -  Im^h  >  0.  By  the  Gronwall  inequality  we  obtain 


W<M4  r||£f-1(r)-£f-2(r)||2dr, 

Jo 


(4.87) 


with  //.|  =  /i3eM2tl,  so 


I^J)l|2+ll«f(0H2  <  M2  /*M4  /r||cf-1(,)-£)V-2(S)||2^T 

JO  Jo 

+  h3  r||£f-1(r)-£f-2(r)||2dr 
4o 


<  (^2/^4 1\  + 


|£rV)^f-2(T)||2dT. 


(4.88) 
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Similarly, 


— Ci<  J*e-^-%v(u»(s))  -  9v(u^mds,m) 

+  ( e  Cl^  S'(<?m(e]V  1(s))  —  dini^l  2{S)))d$i;  (t  )) . 

Using  similar  techniques  to  those  above,  we  obtain  that  for  some  constants  //s ,  inde¬ 

pendent  of  N 


^)ir<^newir+/i6  /*  11^)11^+^  2dr.  (4.89) 

Jo  Jo 


Substituting  (4.87)  into  (4.89)  we  hnd 


\^m2  <  »rmr)-4-\T)\r*, 

Jo 


which  together  with  (4.88)  yields 


(4.90) 


h?m2+\\u*m2+\\i»m2  <  t /tiief-lw-®f-2wir*» 

Jo 

~  f'^'x  '  O!  II'. 2(0, (4.91) 

where  ft  =  /^/Mi  +  !?>  +  ^8-  This  guarantees  that  {uN  (f)},  {u^(t)},  {eA  (0}  are  Cauchy  se¬ 
quences.  Using  the  strong  convergence  of  these  sequences  we  can  take  a  limit  in  (4.77)  and 
(4.80)  to  obtain  the  existence  of  a  weak  solution  of  (4.77)-(4.80)  on  the  interval  [0,U]-  Unique¬ 
ness  of  the  weak  solution  can  be  derived  in  the  usual  way.  Now  the  weak  solution  can  be 
extended  to  the  intervals  [U,  ti+i]  as  in  [6].  Thus  we  proved  that  (4.77)-(4.80)  has  a  unique 
global  weak  solution. 

Remark  4.1  It  is  possible  to  establish  the  global  existence  of  a  weak  solution  under  local  Lip- 
schitz  properties  and  growth  conditions  on  the  nonlinear  functions  ge,gvi,gvd,gin.  Taking  uf 
as  a  test  function  in  (f-II)  and  using  standard  inequalities  we  can  show  that  the  iterates 
(f)},  {eA  (f )}  art  bounded  by  a  constant,  independent  of  N.  Thus  the  computations  (f.81)- 
(4-91)  can  be  repeated  in  this  ball  using  the  local  Lipschitz  property  of  the  nonlinear  functions. 
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